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In this paper, we study an optimal execution problem in the case of uncertainty in 
market impact to derive a more realistic market model. Our model is a generalized 
version of that in [6], where a model of optimal execution with deterministic mar- 
ket impact was formulated. First, we construct a discrete-time model as a value 
function of an optimal execution problem. We express the market impact function 
as a product of a deterministic part (an increasing function with respect to the 
trader's execution volume) and a noise part (a positive random variable). Then, 
we derive a continuous-time model as a limit of a discrete-time value function. We 
find that the continuous-time value function is characterized by an optimal con- 
trol problem with a Levy process and investigate some of its properties, which are 
mathematical generalizations of the results in [6] . We also consider a typical exam- 
ple of the execution problem for a risk- neutral trader under log-linear /quadratic 
market impact with Gamma-distributed noise. 

Keywords: market liquidity, optimal execution, uncertain market impact, Levy 
process, viscosity solution, Hamilton-Jacobi-Bellman equation 



1 Introduction 

The optimal portfolio management problem is central in mathematical finance theory. 
There are various studies on this problem, and recently more realistic problems, such as 
liquidity problems, have attracted considerable attention. In this paper, we focus on market 
impact (MI), which is the effect of the investment behavior of traders on security prices. MI 
plays an important role in portfolio theory, and is also significant when we consider the case 
of an optimal execution problem, where a trader has a certain amount of security holdings 
(shares of a security held) and attempts to liquidate them before the time horizon. The 
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optimal execution problem with MI has been studied in several papers (pQ, [2J, [3J, [I], PS] 
and references therein,) and in [6J such a problem is formulated mathematically. 

It is often assumed that the MI function is deterministic. This assumption means that 
we can obtain information about MI in advance. However, in a real market it is difficult to 
estimate the effects of MI. Moreover, it often happens that a high concentration of unexpected 
orders will result in overffuctuation of the price. The Flash Crash in the United States stock 
market is a notable precedent of unusual thinning liquidity: On May 6th, 2010, the Dow 
Jones Industrial Average plunged by about 9%, only to recover the losses within minutes. 
Considering the uncertainty in MI, it is thus more realistic and meaningful to construct a 
mathematical model of random MI. Moazeni et al. [9] studied the uncertainty in MI caused by 
other institutions by compound Poisson processes, and then studied an optimization problem 
of expected proceeds of execution in a discrete-time setting. They considered the uncertainty 
in arrival times of large trades from other institutions; however, MI functions of decision 
makers themselves were given as deterministic linear functions so that the decision makers 
knew how their own execution affected the market price of the security (the coefficients of 
MI functions were regarded as "expected price depressions caused by trading assets at a unit 
rate"). 

In this paper, we generalize the framework in [6J, particularly considering a random MI 
function. We follow the approach in [6]: we construct a discrete-time model and take a limit 
to derive a continuous-time model and a corresponding value function. We assume that noise 
in the MI function in a discrete-time model is independent of both time and trading volume. 
The randomness of MI in the continuous-time model is described as a jump of a Levy process. 
We study some properties of the continuous-time value function which are mathematical gen- 
eralizations of the results in [6] . In particular, we find that the value function is characterized 
as a viscosity solution of the corresponding Hamilton- Jacobi-Bellman equation (abbreviated 
as HJB) when the MI function is sufficiently strong. We also perform a comparison with the 
case of deterministic MI and show that noise in MI makes a risk-neutral trader underestimate 
the MI cost. This means that a trader attempting to minimize the expected liquidation cost 
is not particularly sensitive to the uncertainty in MI. Moreover, we present generalizations 
of examples in [6J and investigate the effects of noise in MI on the optimal strategy of a 
trader by numerical experiments. We consider a risk-neutral trader execution problem with a 
log-linear /quadratic MI function whose noise is given by the Gamma distribution. 

The rest of this paper is organized as follows. In Section [2J we present the mathematical 
formulation of our model. We set a discrete-time model of an optimal execution problem as our 
basic model and define the corresponding value function. In Section [3j we present our results, 
showing that the continuous-time value function is derived as a limit of the discrete-time one. 
Moreover, we investigate the continuity of the derived value function. Note that the results in 
this section are of the same form as those in [5J. In Section HI we study the characterization 
of the value function as a viscosity solution of the corresponding HJB equation as a direct 
consequence of the result in [6J . In Section |5j we consider the case where the trader must 
sell all the shares of the security, which is referred to as a "sell-off condition." We also study 
the optimization problem under the sell-off condition and show that the results in Section 4 
in [BJ also hold in our model. Section [BJ treats the comparison between deterministic MI and 
random MI in a risk-neutral framework. In Section [7] we present some examples based on the 
proposed model. We conclude this paper in Section [HJ 
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2 The Model 



In this section, we present the details of the proposed model. Let (Q, J 7 , P) be a complete 
probability space. T > denotes a time horizon, and we assume T = 1 for brevity. We 
assume that the market consists of one risk- free asset (cash) and one risky asset (a security). 
The price of cash is always 1, which means that a risk- free rate is zero. The price of the 
security fluctuates according to a certain stochastic flow, and is influenced by sales performed 
by traders. 

First, we consider a discrete-time model with a time interval 1/n. We consider a single 
trader who has an endowment of $ > shares of a security. This trader liquidates the shares 
$o over a time interval [0, 1] considering the effects of MI with noise. We assume that the 
trader sells shares at only times 0, 1/n, . . . , (n — l)/n for n G N = {1, 2, 3, . . .}. 

For I = 0, ... ,n, we denote by S™ the price of the security at time l/n, and we also denote 
Xf = logS*™. Let so > be an initial price (i.e., S$ = so) and Xq = logso- If the trader 
sells an amount at time l/n, the log price changes to X? — g™(i/}™), and by this execution 
(selling) the trader obtains an amount of cash ip^S] 1 exp(— g™(^™)) as proceeds. Here, the 
random function 

^ n (V>, u) = J/{u)g n {il)), V e [0, $o], w G 

denotes MI with noise, which is given by the product of a positive random variable c" and a 
deterministic function g n : [0, $o] — > [0, oo). The function g n is assumed to be non-decreasing, 
continuously differentiable and satisfying g n (0) = 0. Moreover, we assume that (c")j is i.i.d., 
and therefore noise in MI is time- homogeneous. Note that if c" is a constant (i.e., cf = c for 
some c > 0,) then this setting is the same as in [B]. 

After trading at time l/n, X™ +1 and S'[! l _ 1 are given by 

X? + i = y(—; -,X? - g?(W)), S? +1 = e x ^, (2.1) 
V n n / 

where F (t; r, x) is the solution of the following stochastic differential equation (SDE) on the 
filtered space (Q, T , {Tf\, P) 

( dY{t; r, x) = a(Y(t; r, xj)dB t + b(Y(t; r, x))dt, t > r, 
\ Y(r; r, x) = x, 

where (-Bt)o<*<i is standard one-dimensional Brownian motion which is independent of (c™);, 
(J~f)t is its Brownian filtration, and b, a : R — > M are Borel functions. We assume that b 
and er are bounded and Lipschitz continuous. Then, for each r > and there exists a 

unique solution. 

At the end of the time interval [0,1], the trader has an amount of cash W™ and an amount 
of the security where 

W%. 1 = Wr + WS?e-*W), = - w (2.2) 

for I = 0, ... ,n — 1 and Wq =0, ip^ = $ - We say that an execution strategy (^; n )J!T is 
admissible if (ip?)i e -4™(*o) holds, where Al(ip) is the set of strategies (V'Df^o 1 sucn that 
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is J 7 ™ = a {(Bt)t<i/ni c o i ' ' ' i c™_ ^-measurable, ip™ > for each I — 0, . . . , k — 1 and ift™ < <p 

1=0 

almost surely. 

Then, the investor's problem is to choose an admissible strategy to maximize the expected 
utility E[u(W™, <p™, S™)], where u G C is the utility function employed by the investor and C is 
the set of non-decreasing, non-negative and continuous functions on D = K x [0, $o] x [0, oo) 
such that 

u(w,<p,s)<C u (l + \w\ m "+s m *), (w,cp,s)eD (2.3) 

for some constants C u , m u > 0. 

For k = l,...,n, (w,tp,s) G -D and u G C, we define the (discrete-time) value function 

Vk(w,<p,s;u) by 

y>, «; u) = sup E[«(WJ, pj, S£)] 

subject to d2IJ and (Q for Z = 0, . . . , k - 1 and (W™, S£) = (w, </?, s) (for s = 0, we 
set 5" = 0). We denote such a triplet of processes (W™, <p™, )f =0 by S£(u;, y?, s; (^ n )/)> an d 
denote V n (w, 99, s; w) = u(w,<p,s). Then, this problem is equivalent to V™(0, $0, so; u). We 
consider the limit of the value function V£(w, <p, s; u) as n — > 00. 

Let /i : [0, 00) — > [0, 00) be a non- decreasing continuous function. We introduce the 
following condition for g n (ip). 



[A] lim sup 

n ^°°Ve[o,*o] 



0. 



Moreover, we assume the following conditions for (c™)/. 

[Bl] Define 7 n = essinf c™(u;). For any n G N, it holds that 7 n > 0. In addition, 

K x hn) 



n 

holds for x > 0. 



0, n -> 00 (2.4) 



c n + . . . + c n 

[B2l Let [i n be the distribution of ^— . Then, u„ has a weak limit u as n — > 00. 

n 

[B3] There is a sequence of infinitely divisible distributions (p n )n on M such that fi n = fi *p n 
and either 



[B3-a] / x 2 p n (dx) = 0(l/n 3 ) as n -)> 00 



or 

[B3-b] There is a sequence (if n ) n C (0, 00) such that _/\~ n = 0(l/n), p n ((—oo,—K n )) 

= (or p n ((K n , 00)) = 0) and / xp n (dx) = 0(l/n) as n -> 00, where O denotes 

ii 

the order notation (Landau's symbol). 
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Let us discuss condition [Bl]. First, note that 7„ is independent of /, because cf, I = 
0, 1, 2, . . . are identically distributed. Next, we study when the convergence (12. 4p holds. Since 
h is non-decreasing, we see that 

hjxhn) < MOO) 

n n 

where /i(oo) = lim h(Q G [0, oo] (which is well-defined by virtue of the monotonicity of h). 

This inequality tells us that (12 ,4p is always fulfilled whenever h(oo) < oo. In the case of 
h(oo) = oo, we have the following example: 

h(Q = «C P , In = (P,S>0,S< 1/p). (2.5) 

We can actually confirm ( 12.41) by observing 

h(x/j n ) ax p 

— ^ — - = — j- — > 0, n -> oo. 

n n v ° 

Here, the condition [B2] holds only when lim sup 7„ < oo. Indeed, under [B2], we easily see 

n 

that the support of the distribution /x is included in the interval [lim sup j n , oo). Note that 7 n 

n 

of (I2.5p satifies lim sup 7„ < 1 because of the relation 5 < 1/p. 

n 

Since \i is an infinitely divisible distribution, there is some Levy process (subordinator) 
(L t )o<t<i on a certain probability space such that L\ is distributed by /i. Without loss of 
generality, we may assume that (L t )t and (B t ) t are defined on the same filtered space. Since 
(c] 1 ); is independent of (B t )t, we may also assume that (L t )t is independent of (B t )t- Let v be 
the Levy measure of (L t ) t . We assume the following moment condition for v. 

[C] f (z + z 2 )u(dz) < oo. 

i(0,oo) 

Now, we define the function which gives the limit of the discrete-time value function. For 
t e [0, 1] and ip G [0, $o] we denote by At((p) the set of (Jv)o<r<r adapted and caglad processes 
(i.e., left- continuous and having a right limit at each point) (( r )o<r<t such that ( r > for each 
-t 



r G [0, t], / C, r dr < if almost surely and 
Jo 



sup Cr( w ) < oo, (2.6) 

(r,u)e[o,t]xn 

where J> = a{B v , L v ;v < r} V {Null Sets}. Here, the supremum in (12.61) is taken over all 
values in [0,t] x Q. Note that we may replace sup in (I2.6P with esssup. 
For t G [0, 1], (w, <p,s) G D and u G C, we define V t (w, <p, s; u) by 

V t (w,(p,s;u)= sup E[u(W t ,<p t ,S t )} (2.7) 

(Cr)r£A(^) 
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subject to 

dW r = ( r S r dr, 
d<fr = —( r dr, 

dX r = a(X r )dB r + b(X r )dr — g(( r )dL r , S r = exp(X r ) 

and (Wo,ipo,So) = (w,<p,s), where a(s) = sa (log s),b(s) = s|fe(logs) + -cr(logs) 2 | for 

8 > (a(0) = 6(0) = 0) and g(() = [ h((')d('. We denote such a triplet of pro- 

cesses (W r , ip r , S r ) < r < t by E t (w, ip, s; (Cr)r)- Note that Vo(tf , v?, s; tt) = u(w,ip,s). We call 
Vt(w, y?, s; ix) a continuous-time value function. Also note that Vt(iu, s; it) < oo for each 
t G [0, 1] and (w,tp,s) G £>. 



3 Properties of Value Functions 

First, we give the convergence theorem for value functions. 
Theorem 1. For each (w, (p,s) G D, t G [0, 1] and u £ C it holds that 

lim V\lt]( w ^iS,u) = V t (w^,s]u), 

n— »oo 1 J 

where [nt] is the greatest integer < nt. 

According to this theorem, a discrete-time value function converges to Vt(w, <p, s; u) by 
shortening the time intervals of execution. This implies that we can regard Vt(w, tp, s; u) as 
the value function of the continuous-time model of an optimal execution problem with random 
MI. This result is almost the same as in [B], with the exception that the term of MI is given 
as an increment g(( r )dL r . Let 

L t = jt + [ [ zN(dr, dz) 

JO </(0,oo) 

be the Levy decomposition of [L t )t, where 7 > and N(-, ■) is a Poisson random measure. 
Then, g(( r )dL r can be divided into two terms as follows: 

g((r-)dL r = ig(Qdr + flf(Cr-) / zN(dr, dz). 

J (0,00) 

The last term on the right-hand side indicates the effect of noise in MI. This means that noise 
in MI appears as a jump of a Levy process. Using the above representation and Ito's formula, 
we see that when s > 0, the process (S r ) r satisfies 

dS r = &(3 r )dB r + b(S r )dr - \ ig (( r )S r dr + S r . f (1 - e - 9 ^~ )z )N(dz, dr) 

I i(0,oo) 

Regarding the continuity of the continuous-time value function, we have the following 
theorem. 
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Theorem 2. Let u G C. 

(i) If h(oo) = oo, then Vt(w, p, s; u) is continuous in (t, w, tp, s) E [0, 1] x D. 

(ii) If h(oo) < oo, then V t {w, p, s; u) is continuous in (t, w, <p, s) G (0, l]x D and V t {w, tp, s; u) 
converges to Ju(w,tp, s) uniformly on any compact subset of D as t \. 0, where Ju(w,tp, s) is 
given as 



This is also quite similar to the result in [5J, whereby continuities in w,<p and s of the 
value function are always guaranteed, but continuity in t at the origin depends on the state 
of the function h at infinity. When h(oo) = oo, MI for large sales is sufficiently strong (g(() 
diverges rapidly with ( — > oo) to make a trader avoid instant liquidation: an optimal policy is 
"no trading" in infinitesimal time, and thus Vt converges to u as t J, 0. When h{oo) < oo, the 
value function is not always continuous at t = and has the right limit Ju(w,p, s). In this 
case, MI for large sales is not particularly strong (g(() still diverges, but the divergence speed 
is low) and there is room for liquidation within infinitesimal time. The function Ju(w, tp, s) 
corresponds to the utility of the liquidation of the trader, who sells part of the shares of a 
security i[) by dividing it infinitely within an infinitely short time (sufficiently short that the 
fluctuation in the price of the security can be ignored) and obtains an amount p — tp, that is, 



</(0,oo) 

the result. Note that if 7 = and h(oo) < 00, then the effect of MI disappears in Ju(w, <p, s). 
This situation may occur even if E[c£] > Sq (or E[Xi] > £0) for some Sq > 0. 
For the proofs of Theorems HH21 please refer to [5] . 

4 Characterization of the Value Function as a Viscosity 
Solution of a Corresponding HJB 

As presented in [B], we can characterize our value function as a viscosity solution of a 
corresponding HJB when h(oo) = 00. First, we present the Bellman principle (dynamic 
programming principle). Let us define Q t : C — )■ C by Q t u(w, p, s) = V t (w, p, s; u). Then, we 
can easily confirm that Qt is well-defined as a nonlinear operator. The same proof as that for 
Theorem 3 in [6J gives the following proposition. 

Proposition 1. For each r,t G [0, 1] with t + r < 1, (w, <p,s) G D and u G C, it holds that 
Q t+r u(w, p, s) = Q t Q r u(w, y?, s). 

Note that this proposition is also required to prove Theorem El At this stage, we introduce 
the HJB corresponding to our value function. We start by formal calculation to derive the 
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HJB and subsequently show that we can apply Theorems 4-5 in [B] directly to our case, in 
spite of the fact that our value function is a generalization of that in [6] . 

We begin by fixing t G (0, 1] and letting h > be such that t — h > 0. Then, Proposition 
[H leads us to 

V t (w, (f, s; u) = Q h Q t _ h u(w, <p, s; u) = sup E[V t - h (W h , <p h , S h ; u)]. 

(Cr)reA h (ip) 

If we restrict the admissible strategies to constant ones and assume the smoothness of V t (w, if, s; u), 
we obtain 

= sup E[V t _ h (W h , ip h ,S h ;u)} - V t (w,ip, s;u) 

fh 



supE[^ (-^ + Se^j V r (W r , <p r , S r ; u)dr] 



to arrive at 

1 

h 



supE[| ("^ +^ C ) V r (W r ,<p r ,S r ;u)dr] = (4.1) 



by virtue of Ito's formula, where 

Id 2 d 
^v(t,w,tp,s) = -a(s) 2 —v(t,w,Lp,s) + b(s)—v(t,w,Lp,s) 

+C ( s J"^*' w > s ) ~ §^ v ^ t} W} v? ' ~~ 9(C)s-^v(t, w, if, s), 
9(0 = 19(0+ [ (l-e-^>)v(dz). 

J(P,oo) 

We remark that Jzf* is of exactly the same form as in [6]. Letting h — > in ( 14. ip . we obtain 

Q 

— V r (W r , ip r , S r ; u) — sup Jf*V r (W r , ip r , S r ; u) = 0. 

Note that this calculation is nothing but intuitive and formal. However, we can justify the 
following theorem. 

Theorem 3. Assume that h is strictly increasing and h(oo) = oo. Moreover, assume 

hm inf Vt(^s + E;u)-V t (w,<p,s- t u) > Q 
eio e 

for any t G (0, 1] and (to, <p, s) G U. Then, V t (w, ip, s; u) is a viscosity solution of 

d 

—v(t,w,ip,s) — su.p^f < 'v(t,w,(p,s) = 0, (t, w, ip, s) G (0, 1] x U, (4.3) 
ot c>o 

where U = D\dD and D = Kx [0, oo) x [0, oo). 



Proof. We consider the strategy-restricted version of the value function 

V t L (w,p, s;u) = sup E[u(W t ,<fit,St)] 

for L > 0, where A^{tp) = {(( r ) o<r<t £ ^(v 9 ) j su PlCr(^)| < £}■ Then, the standard 
argument (see Section 5.4 in [TT]) suggests that V t L (w, p, s; u) is a viscosity solution of 
d 

—v(t,w,tp,s) — sup JZ**v(t, w, tp, s) = 0, (t, u>, s) £ (0, 1] x [/ 

W 0<C<i 

because of the compactness of the control region [0,L]. Here, we can define h(Q as the 
derivative of g by 

Ho = (V / ^-^1/(^)^(0- 

V -/(0,oo) / 

Obviously, it holds that h(oo) = oo. Hence, we can apply the same argument as in Section 7.6 
in [6] to obtain the assertion by letting L — > oo and using the stability arguments for viscosity 
solutions. ■ 

Since the condition lim(h(()/() > implies lim (/i(C)/C) > 0, we can apply Theorem 5 
in [6] to arrive at the following uniqueness theorem. 

Theorem 4. Assume that a and b are both Lipschitz continuous. Assume the hypotheses in 
Theorem^ and that liminf(/i(£)/£) > 0. If a polynomial growth function v : [0, 1] x D — > R 

is a viscosity solution of (14. 3p and satisfies the boundary conditions 

v(0,w, tp, s) — u(w, tp, s), (w,tp,s)et), 
v(t,w,0,s) = E[u(w,0,Z(t;0,s))], (t,w, s) e [0, 1] x R x [0, 00), (4.4) 
v(t, w, tp, 0) = u(w, tp, 0), (t, w, tp) e [0, 1] x R x [0, 00), 

then V t {w, p, s; u) = v(t, w, tp, s), where 

Z(t;r,s) =exp(F(t;r,logs)) (s > 0), (s = 0). (4.5) 



5 Sell-Off Condition 

In this section, we consider the optimal execution problem under the "sell-off condition" , 
which was introduced in [6]. A trader has a certain quantity of shares of a security at the 
initial time, and must liquidate all of them by the time horizon. Then, the spaces of admissible 
strategies are reduced to 

m^Aiip)-, 5>r 

z=o 

A S °(^) = |(Cr)r £ M<P) I J Crdr 
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Now, we define value functions with the sell-off condition by 

V? 30 {w,<p,s;U) = sup E[U(W£)], 

wr)i6^' so (^) 

V*°{w,<p,s;U) = sup E[U(W t )} 

for a continuous, non- decreasing and polynomial growth function U : K — > R, where 
^, ^ n )f =0 = H£(w, ip, s; (^P)j) and (W r , y> r , S r ) <r<t = H^w, s; (&.),.). 
The following theorem is analogous to Theorem 7 in [6]. 

Theorem 5. V^°(w, ip, s; U) = V t (w, ip, s; u), where u(w, <p, s) = U(w). 

Proof. The relation V^°(w, ip, s; U) < Vt(w, tp, s;u) is trivial, so we show only the assertion 
V t (w, <p, s; U) > V t (w, ip, s; u). Take any (( r ) r G A t (tp) and let (W r , <p r , S r ) r = Hi(iy, ip, s; (Cr)r)- 
Also, take any 5 G (0,£). We define an execution strategy (C*)r £ ^^(v 9 ) by Cr = Cr ( r G 
[0, *-<$]), p*-«/<* (r G (t-M)- Let (W r 5 , ^, S r 5 ) r = S*(u;, <p, s; (Or)- Then, we have 
W t _ a = W/^ < W/, arriving at E[*7(Wi_i)] < E[U(Wf)] < V t so (w,<p,s;U). Letting 5 | 0, 
we obtain E[C/(Wt)] < V^°(w,(p, s;U) by using the monotone convergence theorem. Since 
(Cr)r £ A t ((p) is arbitrary, we obtain the assertion. ■ 

By Theorem we see that the sell-off condition / ( r dr = <p does not introduce changes 

Jo 

in the (value of the) value function in a continuous-time model. Thus, although the value 
function in a discrete-time model may depend on whether the sell-off condition is imposed, in 
the continuous-time model this condition is irrelevant. 

The following is also similar to Theorem 7 in [B], which is a version of Theorem [T] with the 
sell-off condition. 

Theorem 6. For any (w, <p, s) G D, 

hm V, n J°(w, <p, s- U) = V t SO {w, <p, s; U) (= V t (w, cp, s; U).) 

Proof. We may assume t > 0. Take any 5 G (0,t) and let n > 1/8. Then, each strategy 
in A[ n ( t _ s ^(<p) can always be extended to the one in A?^ (<p) by liquidating all remaining 
inventory in the last period. Thus, we see that for n > 1/8 

V[n(t-5)]( w , <P, s ; u ) < ^°(w> <P, s; U) < Vfc t] (w, <p, s; u). (5.1) 

By Theorem [TJ we obtain 

lim V l nf t -s)](w,(p t 8)u) = Vt-s(w,<p,8-,u), lim V£*(w, <p, s; u) = V t (w : <p } s; u). (5.2) 

By (15. ip . (15. 2p . and Theorem [2J we obtain the assertion. ■ 

Analogously to Theorem 8 in [6J, a result similar to Theorem 3 in [8] holds when g(Q is 
linear: 
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Theorem 7. Assume g(() = a>o( for a > 0. 

(i) Vf°(u>, s; Efi = V? (w + s, e~ m s; U ) , where 

V 7"o 



= sup E[U(W t )] 

(Tp r ) r eA t {tp) 



s.t. dS r = e- m 4(S r e m ^dr + e-™^a{S r e ia °^)dB r - S r -dG r , 

(iVT r = dS r , 

So = s, W = w 



and 



where 



A t (<p) = \(<P- I Cvdv) ; (Cr)o<r<t e A S °(^) 



G r = / (1 -e- ao( °- z )N{ds,dz). 

JO </(0,oo) 



(ii) If U is concave and b(s) < for s > 0, then 



V f bU (w,(f,s;U) = U [w+ s . (5.3) 

V 7"o / 

Proof. We can easily confirm assertion (i) by applying Ito's formula to S r and lV r . By a 
similar argument to that in Section 7.9 in [6], we obtain 



E[U(W t )} < U (w + jf E 



■6(,S r e w 0- / S r (l - e~ ao ^- z )v{dz) 



7«o 7(0,(30) 7«o 



for any (<£> r ) r G -4.t(y?) by virtue of the Jensen inequality. Since b is non-positive, the function 
U is non-decreasing, and the terms 

J _ e ~7aoVr g70W r _ I I _ e ~a ( r -z 

are all non-negative, we see that E[£^(W t )] < U(w) for any (jp r ) r e At(<f), which implies 
Vf(w,s) < U(w). The opposite inequality V^(w,s) > U(w) is obtained, similarly to the 
result in Section 7.9 in [6] and the proof of Proposition 11 in [5]. ■ 

We note that the assertion (ii) is the same as Theorem 3 in [5] , and in this case we can obtain 
the explicit form of the value function. The right-hand side of ( 15. 3 j) is equal to Ju(w, <p, s) for 
u(w, ip,s) = U (w) and the nearly optimal strategy for V t ip, s; U) = V t (w, ip, s; u) is given 
by (13. ip . This implies that whenever considering a linear MI function, a risk-averse (or risk- 
neutral) trader's optimal strategy of liquidating a security that has a negative risk-adjusted 
drift is nearly the same as block liquidation (i.e. selling all shares at once) at the initial time. 
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6 Effect of Uncertainty in MI in the Risk-Neutral Frame- 
work 



The purpose of this section is to investigate how the noise in the MI function affects the 
trader. Particularly, we focus on the case where the trader is risk-neutral, i.e. u(w, tp, s) = 
u w (w,tp, s) = w. 

First, we prepare a value function of the execution problem with a deterministic MI function 
to perform a comparison with the case of random MI. Let gJXVO be as in Section [3j and set 
g n (ijj) = E[fl , J l (V')] = E[ c ™] 9n(i>)- Note that when E[c"] = 1, the function g n is equivalent to g n . 
We denote the discrete-time value function by V£{w, tp, s; u) with an MI function g n . Theorem 
1 in [6j implies that under [A] the function V^Aw, tp, s; u) converges to the continuous-time 
value function V t (w, tp, s; u), which becomes the same as in (12. 7p by replacing g(() and L t with 
7<?(C) and t, i.e., the SDE for (X r ) r is given as 

dX r = a(X r )dB r + b(X r )dr - jg(( r )dr, 

where 7 = E[Ia]. 

Theorem 8. V£(w, tp, s; u KN ) > V£(w, tp, s; n RN ). 

Proof. Take any G A^(<p) and let (W, n ,^,S^)j = E%(w, tp, s; (^ re )0 be the triplet for 

V^(w, tp, s; Wrn)- Then, the Jensen inequality implies 

fe-i 

eto = w+"£nws?eM-9nm)} 
1=0 
fe-i 

= w + m i S? exp(- E[$mgn(1>?))] 

1=0 
k-1 

< w + ^E[WS?E[exp{-<?9 n (W))\?r]] < V k n (w,tp,s;u RN ). 
1=0 

Since (ipi')i is arbitrary, we obtain the assertion. ■ 

The above proposition immediately leads us to 

V t (w, tp, s; m rn ) > V t (w, tp, s; w RN ). (6.1) 

This inequality shows that noise in MI is welcome since it decreases the liquidation cost for a 
risk-neutral trader. 

For instance, we consider the situation where the trader estimates the MI function from 
historical data and tries to minimize the expected liquidation cost. Then, a higher sensitivity 
of the trader to the volatility risk of MI results in a lower estimate for the expected proceeds of 
the liquidation. This implies that taking into consideration the uncertainty in MI makes the 
trader prone to underestimating the liquidation cost. Thus, as long as the trader's target is the 
expected cost, considering the uncertainty in MI is not an incentive for being conservative with 
respect to the unpredictable liquidity risk. In Section [TJ we present the results of numerical 
experiments conducted in order to simulate above phenomenon. 
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7 Examples 



In this section, we show two examples of our model, both of which are generalizations of 
the ones in [6]. 

Motivated by the Black-Scholes type market model, we take b(x) = —fi and cr(x) = a for 
some constants \i, a > and assume p, = fj, — cr 2 /2 > 0. We also assume that a trader has a 
risk-neutral utility function u(w,tp,s) = urn(u') = if. In this case, if there is no MI, then a 
risk-neutral trader is afraid of decreasing the expected stock price, and hastens to liquidate 
the shares. 

We consider MI functions which are (log-) linear and (log-) quadratic with respect to liqui- 
dation speed, and assume noise distributed by the Gamma distribution. 
For the noise part of MI, we set 7 n = 7 for n6N and 

P(c™ — 7 G dx) = Gamma(ai/n, nf3x){dx) 

= 7T7 , / g ttl/w ~ 1 e" a/(w/9l) l(o 00 )(g)da;, 

T{a 1 /n){n(5 1 ) a ^ n 1 ' n ' 

where T{x) is the Gamma function. Here, ati, flx, and 7 > are constants. 

For the deterministic part of MI, we consider two patterns such that g n {ip) = &oip and 
<7n(V0 = naoip 2 for ct > 0. In each case, assumptions [A], [B1]-[B3] and [C] are satisfied. The 
corresponding Levy measure is 

V {d z ) = ^ e - z '^l M {z)dz. 



7.1 Log-Linear Impact &; Gamma Distribution 

Theorem [7] directly implies the following theorem. 
Theorem 9. It holds that 

I _ e -7tw 

V t (w,<p,s;u VN ) = w-\ s (7.1) 

7«o 

for each t G (0, 1] and (w, if, s) G D. 

The implication of this result is the same as in [6]: the right-hand side of (17. ip is equal to 
Ju(w, ip, s) and converges to w + <ps as a I or 7 | 0, which is the profit gained by choosing 
the execution strategy of the so-called block liquidation, where a trader sells all shares ip at 
t — when there is no MI. Therefore, the optimal strategy in this case is to liquidate all shares 
by dividing infinitely within an infinitely short time at t — (we refer to such a strategy as a 
nearly block liquidation at the initial time). Note that the jump part of MI 

g((r-) [ zN(dr,dz) 
J (0,00) 

does not influence the value of V t (w, ip, s; mrn)- 
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7.2 Log-Quadratic Impact &; Gamma Distribution 

In [6] , we obtained a partial analytical solution to the problem: when tp is sufficiently small 
or large, we obtain the explicit form of optimal strategies. However, noise in MI complicates 
the problem, and deriving the explicit solution is more difficult. Thus, we rely on numerical 
experiments. Owing to the assumption that the trader is risk-neutral, we may assume that 
an optimal strategy is deterministic. Here, we introduce the following additional condition. 

P] 7 > 

In fact, we can replace our optimization problem with the deterministic control problem 
fit, ip) = sup/ exp I- q(Cv)dv) Crdr, 

(Cr)r J0 \ J0 J 

?(C) = P> + 7"oC 2 + "i log(« /3iC 2 + 1) 

for a deterministic process (( r ) r under the above assumption. We quote Theorem 5 in [5] 
below. 

Theorem 10. V t (w, <p, s; mrn) = w + sf(t, <p) under [D] . 

Here, note that Theorems [3HH imply the following characterization. 
Proposition 2. f(t,tp) is the viscosity solution of 

9 -f ~ sup (/ + C ( 1 - -0-f) - 7«oC 2 / + «i log(« /3iC 2 + 1)1=0 (7.2) 



dt ? > i V dip' 
with the boundary condition 

f{0,<p) = <p, /(*,0) = 0. (7.3) 

Moreover, if f is a viscosity solution of ( 17. 2p — f !7.3j) and /ias a polynomial growth rate, then 
/ = /• 

It is difficult to obtain an explicit form of the solution of (l7.2p -( 17T3"j) . Instead, we solve this 
problem numerically by considering the deterministic control problem (ip) in the discrete- 
time model for a sufficiently large n: 



fc-i , , 

sup exp x ~ ~ Yl { n ^Wf + — log(n 2 a /3i W ; 

WDf^cM*. , =0 \ n m =o L n 

We set each parameter as follows: a = 0.01, t — 1,/j = 0.05, w = 0, s = 1 and n = 500. For 
we examine three patterns for <p = 1, 10, and 100. 
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7.2.1 The case of fixed 7 

In this subsection, we set 7 = 1 to examine the effects of the shape parameter oi\ of the 
noise in MI. Here, we also set j3\ = 2. As seen in the numerical experiment in [6], the forms 
of optimal strategies vary according to the value of ip. Therefore, we summarize our results 
separately for each <p. 

Figure [T] shows the graphs of the optimal strategy (( r )r and its corresponding process 
{<Pt)t of the security holdings in the case of <p = 1, that is, the number of initial shares of 
the security is small. As found in [6], if there is no noise in the MI function (i.e., if a± = 0), 
then the optimal strategy is to sell up the entire amount at the same speed (note that the 
roundness at the corner in the left graph of Figure [T] represents the discretization error and is 
not essential). The same tendency is found in the case of ot\ = 1, but in this case the execution 
time is longer than in the case of ot\ = 0. When we take a\ = 3, the situation undergoes a 
complete change. In this case, the optimal strategy is to increase the execution speed as the 
time horizon approaches. 

When the amount of the security holdings is 10, which is larger than in the case of tp — 1, 
the optimal strategy and the corresponding process of the security holdings are as shown in 
Figure EJ In this case, a trader's optimal strategy is to increase the execution speed as the end 
of the trading time approaches, which is the same as in the case of p = 1 with a\ = 3. Clearly, 
a larger value of ot\ corresponds to a higher speed of execution closer to the time horizon. We 
should add that a trader cannot complete the liquidation when a\ = 3: However, as mentioned 
in Sectional we can choose a nearly optimal strategy from Af°(<p) without changing the value 
of the expected proceeds of liquidation by combining the execution strategy in Figure [2] (with 
«i = 3) and the terminal (nearly) block liquidation. See Section 5.2 of [6] for details. 

When the amount of the security holdings is too large, as in the case of tp = 100, a trader 
cannot complete the liquidation regardless of the value of a±, as Figure [3] implies. This is 
similar to the case of tp — 10 with a± — 3. The remaining amount of shares of the security 
at the time horizon is larger for larger noise in MI. Note that the trader can also sell all 
the shares of the security without decreasing the profit by combining the strategy with the 
terminal (nearly) block liquidation. 

7.2.2 The case of fixed 7 

In the above subsection, we presented a numerical experiment performed to compare the 
effects of the parameter ot\ by fixing 7. Here, we perform numerical comparison from a different 
viewpoint. 

The results in Section [6] imply that taking the uncertainty in MI into account makes a risk- 
neutral trader optimistic about the estimation of liquidity risks. To obtain a deeper insight, 
we investigate the structure of the MI function in more detail. In Theorems [2](ii) and O 
the important parameter is 7, which is the infimum of L\ and is smaller than (or equal to) 
E[£i]- We can interpret this as a characteristic feature whereby the (nearly) block liquidation 
eliminates the effect of positive jumps of {L t )t- However, there is another decomposition of L t 
such that 
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where 

7 = 7 + 

'(0,00) 



7 = 7+/ zv{dz) = E[Li 
J (0,00) 



and 



N(dr, dz) = N(dr, dz) - / zv{dz)dr. 

'(0,00) 



This representation is essential from the viewpoint of martingale theory. Here, N(-, •) is the 
compensator of N(-, ■) and 7 can be regarded as the "expectation" of the noise in MI. Just 
for a risk- neutral world (in which a trader is risk-neutral), as studied in Section [61 we can 
compare our model with the case of deterministic MI functions as in [Hj by setting 7 = 1. 
Based on this, we conduct another numerical experiment with a constant value of 7. 
Note that in our example 



7 = 7 + cuiA (7-4) 



and 

1, 
t 



Var( / / zN(dr, dz) ) = atfl (7.5) 
'(0,00) 



hold. Here, (j7.4p (respectively, (j7.5j) ) corresponds to the mean (respectively, the variance) of 
the noise in the MI function at unit time. Comparisons in this subsection are performed with 
the following assumptions. We set the parameters (3± and 7 to satisfy 

7 + = 1, atiPi = 0-5- 

We examine the cases of ot\ = 0.5 and 1, and compare them with the case of 7 = 1 and ot\ = 0. 

Figure H] shows the case of tp = 1, where the trader has a small amount of security holdings. 
Compared with the case in Section l7.2.1l the forms of all optimal strategies are the same; that 
is, the trader should sell the entire amount at the same speed. The execution times for a.\ > 
are somewhat shorter than for a± = 0. 

Figure [5] corresponds to the case of tp — 10. The forms of the optimal strategies are similar 
to the case of tp — 10, ax — 0, 1 in Section [7.2. II Clearly, the speed of execution near the time 
horizon increases with increasing a.\. 

The results for p> = 100 are shown in Figure El The forms of the optimal strategies are 
similar to the case of tp = 100 in Section 17.2.11 However, in contrast to the results in the 
previous subsection, the remaining amount of shares of the security at the time horizon is 
smaller for larger ot\. 

Finally, we investigate the total MI cost introduced in [7]: 

TC(* ) = -.o g ii%^. (7 . 6) 

When the market is fully liquid and there is no MI, then the total proceeds of liquidating $ 
shares of the security are equal to $ S , because the drift of the price process is negative and the 
expected price decreases as time passes (the optimal strategy of a risk-neutral trader is block 
liquidation at the initial time). On the other hand, in the presence of MI, the (optimal) total 
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proceeds decrease to V T (0, $o, s) = § s x exp(-TC($ )). Thus, the total MI cost TC($ ) 
denotes the loss rate caused by MI in a risk- neutral world. 

Figure [7] shows the total MI costs in the cases of $o = 1 an d 10. Here, we omit the case of 
$o = 100 because the amount of shares of the security is too large to complete the liquidation 
unless otherwise combining terminal block liquidations (which may crash the market). In both 
cases of $ = 1 and 10, we find that the total MI cost decreases by increasing a\. Since the 
expected value 7 of the noise in MI is fixed, an increase in ot\ implies a decrease in 7 and f3\. 
Risk- neutral traders seem to be sensitive toward the parameter 7 rather than toward a±, and 
thus the trader can liquidate the security without concern about the volatility of the noise in 
MI. Therefore, the total MI cost for ol\ > is lower than that for a,\ = 0. 

8 Concluding Remarks 

In this paper, we generalized the framework in [6] and studied an optimal execution problem 
with random MI. We defined the MI function as a product of an i.i.d. positive random variable 
and a deterministic function in a discrete-time model. Furthermore, we derived the continuous- 
time model of an optimization problem as a limit of the discrete-time models, and found out 
that the noise in MI in the continuous-time model can be described as a Levy process. Our 
main results discussed in Sections [3H5] are almost the same as in [6]. 

When considering uncertainty in MI, there are two typical barometers of the "level" of 
MI: 7 and 7. By using the former parameter 7, we can decompose MI into a deterministic 

part jg(( t )dt and a pure jump part g(Ct-) / zN(dt, dz). Then, the pure jump part can be 

</(0,oo) 

regarded as the difference from the deterministic MI case studied in [6]. On the other hand, 
as mentioned in Sections [6] and [TJ the latter parameter 7 is important not only in martingale 
theory but also in a risk-neutral world. Studying 7 also provides some hints about actual 
trading practices. Regardless of whether we accommodate uncertainty into MI, it may result 
in an underestimate of MI for a risk-neutral trader. The risk-neutral setting (i.e. the utility 
function is set as mrn(w) = w) is typical and standard assumption in the study of the execution 
problem, especially in the limit-order-book model (see, for instance, [I]). 

Studying the effects of uncertainty in MI in a risk-averse world is also meaningful. As 
mentioned in Section [SJ when the deterministic part of the MI function is linear, the uncertainty 
in MI does not significantly influence the trader's behavior, even when the trader is risk-averse. 
In future work, we will investigate the case of nonlinear MI. 

Finally, in our settings, the MI function is stationary in time, but in the real market, the 
characteristics of MI change according to the time zone. Therefore, it is meaningful to study 
the case where the MI function is inhomogeneous in time. This is another topic for future 
work. 
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Figure 3: Result for (p = 100 in the case of fixed 7. Left : The optimal strategy ( r . Right : 
The amount of security holdings tp r . 
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Figure 4: Result for ip — 1 in the case of fixed 7. Left : The optimal strategy ( r . Right : The 
amount of security holdings ip r . 




Figure 5: Result for (p — 10 in the case of fixed 7. Left: The optimal strategy ( r . Right : The 
amount of security holdings </? r . 
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Figure 6: Result for ip = 100 in the case of fixed 7. Left : The optimal strategy ( r . Right : 
The amount of security holdings ip r . 
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Figure 7: Total MI cost TC($o) f° r a risk-neutral trader. Left: the case of $o = 1. Right: 
the case of $o = 10. The horizontal axes denote the shape parameter a\ of the Gamma 
distribution. 
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